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We report on a theoretical study of the commensurability oscillations in a quasi-two-dimensional 
electron gas modulated by a unidirectional periodic potential and subjected to tilted magnetic fields 
with a strong in-plane component. As a result of coupling of the in-plane field component and the 
confining potential in the finite-width quantum well, the originally circular cyclotron orbits become 
anisotropic and tilted out of the sample plane. A quasi-classical approach to the theory, that relates 
the magneto-resistance oscillations to the guiding-center drift, is extended to this case. 
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I. INTRODUCTION 

An isotropic quasi-two-dimensional electron system 
with circular Fermi contours undergoes the Lifshitz phase 
transition 1 under the influence of a strong in-plane mag¬ 
netic field. The topology of Fermi contours changes. The 
deviation from the circular shape depends on the mag¬ 
netic field strength and the form of the confining poten¬ 
tial. In systems with two occupied subbands the excited 
subband is emptied at a certain critical in-plane field, 
and the corresponding second Fermi loop disappears. 

The Fermi contours acquire an asymmetric egg-like 
shape in an a sym metric triangular potential at the 
hetero-interface . 2 3 In wide quantum wells and double 
wells with a single occupied subband the Fermi contours 
resemble the Cassini oval: as the in-plane magnetic field 
increases the elongated convex curves acquire the con¬ 
cave peanut-like shape. At high in-plane fields the single 
Fermi line is split into two partsPEl 

The deformation of a Fermi contour shape can be char¬ 
acterized by a single experimentally measurable quantity, 
the magnetic-field-dependent cyclotron massP^Its field- 
dependence can be studied e.g. by the cyclotron reso¬ 
nance in the infrared region of the optical spectra®®! 
or the temperature damping of Shubnikov-de Haas 
oscillations.®®® The closely related magnetic-field de¬ 
pendence of the density of states is reflected in a resis¬ 
tance oscillatio n m easured as a function of the in-plane 
magnetic field lUJlfl 

More details about the size and shape of cyclotron 
orbits can be gained from the magneto-electron focus¬ 
ing experime nts a nd the commensurability oscillations 
measurement ! 16 17 

The commensurability oscillations 18 21 - oscillations 
of the magnetoresistance, measured at low temperatures 
and in a low perpendicular magnetic field in the presence 
of a weak modulation potential - are periodic in the in¬ 
verse field. Their period reflects the commensurability of 
the cyclotron orbit diameter and the modulation period. 

The first usage of the commensurability oscillations 
measurement in strong in-plane magnetic fields was to 


confirm the distortion of the Fermi contour to the egg- 
like shape . 17 Two experimental arrangements were ex¬ 
amined — with a lattice vector of the unidirectional lat¬ 
eral superlatices either parallel or perpendicular to the 
in-plane magnetic-field component. The observed results 
were consistent with the theoretical prediction. 

The influence of the Fermi loop egg-like deformation 
on the chaotic electron dynamics in a two-dimensional 
antidot lattice was studied both theoretically and exper¬ 
imentally. Reasonable agree ment b etween the theory and 
the experiment was achieved . 22 23 

Recently, the commensurability oscillations were inves¬ 
tigated in detail in a wide double hetero-junction well 
with an occupied bonding subband and a unidirectional 
modulation potential .® 11 The caliper dimensions of the in¬ 
plane field distorted Fermi contours obtained from the 
experimental data were compared with the results of 
the first-principle self-consistent calculation. An overall 
semi-quantitative agreement was achieved between the 
experimental and the theoretical results. However, a sys¬ 
tematic discrepancy was found between the observed and 
the calculated elongation of the Fermi contour for the 
case of a lattice vector parallel to the in-plane magnetic 
field. 

To shed light on this apparent discrepancy between 
the theoretical and the experimental findings, we try to 
extend the quasi-classical approach, that relates the mag- 
netoresistance oscillations to the guiding-center drift, to 
the case of cyclotron orbits which are anisotropic and 
tilted out of the sample plane. 


II. CYCLOTRON ORBITS IN TILTED 
MAGNETIC FIELD 

Let us consider a single bonding subband of a symmet¬ 
ric wide well or a double well. 3 !® Assuming the in-plane 
magnetic field parallel to y-axis, the Fermi contour in the 
k x ~ k y plane is described by an expression 
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where Ep denotes the Fermi energy. Only the k x - 
dependence of the energy is influenced by the in-plane 
field, By , while the harmonic dependence of the energy 
on the wave vector component k y remains untouched. 
With increasing B y the curvature of E(k x ) decreases for 
k x close to k x = 0 and, at a certain value of B y , be¬ 
comes negative. A local maximum develops at k x = 0, 
accompanied by two new minima positioned symmetri¬ 
cally around it. The corresponding Fermi contours re¬ 
semble the Cassini ovals: the convex curves elongated 
in fc x -direction acquire the concave peanut-like shape at 
higher fields with the width of a peanut ‘waist’ at k x = 0 
shrinking to zero as B y increases. 

The quasi-classical theory predicts that in a weak 
perpendicular magnetic field B z an electron is driven 
around a Fermi contour by the Lorentz force with the 
velocity components given by v x = l/HdE(k x )/dk x and 
v y = Hky/m. The related real-space cyclotron orbits are 
similar in shape, but scaled by £ 2 = h/{\e\B z ) and ro¬ 
tated by 7t/2, as v x = £ 2 dk y /dt and v y = —£ 2 dk x /dt. 

From these expressions we can calculate the in-plane 
field-dependent period of the cyclotron motion T, or 
equivalently the cyclotron frequency lo c = \e\B z /m c , 
where the cyclotron effective mass m c is related to the 
shape of the Fermi contour by 
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From here we obtain the density of states g: 

m c 

9= ^h 2 ' 
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Eq. © also determines the relation between the electron 
concentration N and the Fermi energy Ep through N = 
gEp. In zero in-plane field this relation reduces to Ep = 
h 2 kp/2m , where the Fermi wave vector kp is given by 
kp = 

Thus, assuming that the coordinates of a guiding cen¬ 
ter in x—y plane are X and Y, the time dependence of an 
electron motion along the real-space orbits is described 
by equations 


x{t) = X + £ 2 z k y (t), y(t) =Y - £ 2 z k x (t). (4) 


Due to the finite well width d the real-space orbit has 
also a z-coordinateP 

z(t) = £ 2 y (k x - , (5) 


where £ 2 = H/(\e\B y ). 

Two new important features appear when compared 
with previously discussed zero-in-plane-field cyclotron or¬ 
bits. 

First, the velocity of an electron moving along an or¬ 
bit is no longer constant and depends on its position on 
the Fermi contour. Most pronounced changes go for the 


so called turning points, i.e., the points on the Fermi 
contour where one of the velocity components drops to 
zero. While in the zero-field case there were only two 
pairs of turning points at crossings of the Fermi contour 
with the k x and k y axes, new turning points appear when 
the Fermi contour becomes concave. Their position on 
a Fermi line is related to field-induced new minima of 


E(k x ). 

Second, as follows from Eqs.(4|5) the cyclotron orbit is 
tilted. The z-coordinate of an electron does not remain 
in the center of a well, and the electron moves from one 
side of a well to the other when traveling along a real 
space path. 


III. COMMENSURABILITY OSCILLATIONS 

The quasi-classical description of the commensurabil- 
ity oscillations - the magneto-resistance oscillations due 
to the presence of a weak unidirectional potential - is 
based on the solution of the Boltzmann equation in the 
relaxation time approximation 21 . The oscillations are pe¬ 
riodic in 1/B Z and their amplitude is proportional to the 
mean square drift velocity averaged over the positions of 
the guiding centers. The drift velocity of a given guiding 
center can be calculated by the time-average of an oscil¬ 
lating electric field seen by an electron traveling around 
the center along the real-space cyclotron orbit. 

In zero in-plane magnetic field the cyclotron orbits are 
circles with a radius R = £ 2 kp. If we assume a sinusoidal 
potential with the lattice vector a oriented along the y- 
axis, the decisive contributions to the drift velocity are 
obtained at the two turning points Y + R and Y — I?, 
where an electron stops for a while, before the velocity 
component v y changes sign. 

It was shown®! that the frequency / of commensurabil- 
ity oscillations is proportional to the radius of the Fermi 
contour kp by 


/ = 



( 6 ) 


i.e., the frequency / measures the caliper dimension of 
the Fermi contour. 

In the following we discuss why the above formula can¬ 
not be simply extended to the case of non-circular Fermi 
contours by replacing kp by the Fermi contour calipers, 
which determine the position of the turning points on the 
real space orbits. 

In finite B yi the time-averaged electric field of the orbit 
with the guiding center coordinate Y is given by 

£{Y) = j £o (z{t))cos (^-y{t)j dt. (7) 


The corresponding drift velocity reads thrift = £(Y)/B Z , 
the dependence of y and z coordinates on time is de¬ 
scribed through the time dependent k x by Eqs.(4|5). 
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When compared with the zero-field case, the main dif¬ 
ference stems from the following facts. 

There are four more turning points for concave Fermi 
contours where the velocity components drop to zero, as 
shown in the Fig. (jT|). 

Moreover, we assume that the source of the unidirec¬ 
tional periodic modulation is only on one side of a quan¬ 
tum well and therefore the modulation potential is partly 
screened by the quasi-two-dimensional electron layer of a 
finite width. Then the strength of the oscillating electric 
field is 2 -dependent and need not be the same on both 
sides of a well. 



FIG. 1. The caliper dimensions of the concave Fermi contour 
at By = 8T. Dots mark the turning points on the related 
cyclotron orbits rotated by 7t/2 and multiplied by t 2 z . 


As mentioned above, the transport coefficients ob¬ 
tained by the solution of the Boltzmann equation, or 
from the linear response theoryP^hl depend on the mean 
square of the drift velocities. It was shown by both meth¬ 
ods that A Q yy , the change of the magneto-resistance in 
the presence of a weak unidirectional potential, is pro¬ 
portional to the mean square drift velocity averaged over 
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Here r is the relaxation time, qq and cto are the resistivity 
and conductivity at B z = 0, respectively. 

Similar expressions can be written for the lattice vector 
oriented along the ^-direction. 


IV. MODEL CALCULATION 

To illustrate the relevance of the above extension of 
the quasi-classical description of commensurability os¬ 
cillations we present the results of a numerical calcula¬ 
tion based on the simple tight-binding model of a double 
weiiJUHH This model captures quite well the essence of 


physics of double wells and wide quantum wells subjected 
to a strong in-plane magnetic field. 

We consider two strictly two-dimensional electron lay¬ 
ers in very narrow quantum wells separated by a bar¬ 
rier. The model is characterized by two parameters, 
the interlayer distance d and the hopping integral th- 
Such a system represents the first step from strictly two- 
dimensional to three-dimensional structures. 

The diagonalization of a 2 x 2 matrix yields the /in¬ 
dependence of the bonding subband energy in the form 


+ 
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where S is the abbreviation for h\e\dB y k x /2m. The Fermi 
energy Ep corresponding to the fixed concentration of 
carriers N is determined by Eqs. ( |2|3[ ). Here we choose 
N = 8 x 10 11 cm -2 , d = 20nm and th = 2meV. A few 
examples of cyclotron orbits calculated based on these 
parameter are presented in the Fig [2] 



FIG. 2. Projections of cyclotron orbits on the x — y and y — z 
planes, B z ~ 0.07 T. a) B y = 0T. The orbit is circular and 
lies in the middle of a double well at 2 = 0. b) B y = 3T. 
The orbit is on the border between the convex and concave 
curve. The 2 -coordinate changes sign at y = 0. c) B y = 8T. 
The concave shape of the orbit is fully developed. An electron 
spends most time near the opposite sides of the structure at 
2 = ±d/2. d) By = 9.75 T. The orbit is close to splitting. 
The change of sign of the 2 -coordinate near y = 0 is very 
sharp. 


The independence of commensurability oscillations is 
given by Eqs.Q and |8]). These expressions were evalu¬ 
ated by numerical integration for the fixed in-plane field 
components B y between 0 and 9.75 T. The domain of 
B z was chosen between 0.05 and 0.35 T. This is con¬ 
sistent with a range of fields commonly used in experi¬ 
ments - below the minimum field the oscillations are usu¬ 
ally damped, above the maximum field the Shubnikov- 
de Hass oscillations dominate. Note that this choice is 
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not an inherent property of the quasi-classical approach. 
But, as it will be shown later, it can slightly influence 
interpretation of our model calculation. 

To test the accuracy of our numerical method, be¬ 
sides of (Wd rift ) we also calculate the averaged drift ve¬ 
locity (udrift)- It defines the macroscopic current e(udrift) 
through the sample which must, of course, disappear in 
the thermodynamic equilibrium. 

We first test our model on two examples: the stan¬ 
dard circular Fermi contour at B y = 0 T and a strongly 
anisotropic Fermi contour at B y = 8T. The output of 
our model calculation is presented in the Fig(3] 

The results for a circular Fermi contour are shown in 
the Fig|3]a). They agree with results published pre¬ 
viously, the oscillations reach zero in each period and 
their frequency satisfies the commensurability condition 
/ = kp. given by Eq. (J gJ) . The period of oscillations is 
the same for a lattice vector a oriented along the x and 
y- axes. The reason is that an electron traveling along 
the circular orbit stays all the time in the middle of the 
structure, z = 0, and, therefore, the z-dependence of a 
modulation potential plays no role. 

The case of anisotropic Fermi contours is quite differ¬ 
ent. The dependence of oscillations on the lattice vector 
orientation is strong and the z-dependence of the modu¬ 
lation potential can play an important role. 



FIG. 3. Commensurability oscillations as a function of B z . 
a) By = 0T. The orbit is circular and the same form of os¬ 
cillations for a || a: and a || y is obtained. The z-dependence 
of a modulation potential plays no role, b) B y = 8T, a || x. 
Again, the z-dependence of a modulation potential plays no 
role, as an electron feels the same electric field at x and —x. 
c) By = 8 T, a || y and a z-independent modulation potential. 
Two frequencies of oscillations are resolved, their aplitude 
reaches zero, d) B y = 8 T, a || y and a z-dependent mod¬ 
ulation potential. An electron feels different electric fields at 
y and —y. The amplitude of oscillations does not drop to zero. 
The weights of two frequencies in the frequency spectrum are 
changed. 


Let us first consider the lattice vector a parallel to 
the axaxis. As in the previous case, B y = 0T, the z- 
dependence of the modulation potential has almost no 
influence on the resulting shape of commensurability os¬ 
cillations, shown in the Fig{3]b). The reason is that the 
z-coordinate depends only on k x and an electron feels the 
same electric field at x and — x. The fast Fourier trans¬ 
form (FFT) of oscillations yields three frequencies related 
to four turning points which correspond to ± k,F, y max, 
two turning points corresponding to ± kp, y min, and their 
average. 

Now we turn our attention to a || y and a z-independent 
modulation potential. The calculated results are pre¬ 
sented in the Fig{3] c). The FFT of oscillations yields 
two frequencies close to fi = pp- kp x and / 2 = rf kp x . 
As mentioned above their deviations from the exact val¬ 
ues and their weights in the frequency spectrum depend 
slightly on the chosen domain of B~. The frequency /i is 
related to the turning points ± shown in FigfTl The 
frequency / 2 appears for the first time at the critical field 
B y for which the Fermi contour becomes concave, in our 
case at B y = 3 T. The weight of / 2 in the frequency spec¬ 
trum increases with increasing B y and reaches maximum 
at a critical field B y = 9.75 T above which the Fermi con¬ 
tour splits. The reason is that at this field both velocity 
components v x and v y drop to zero at k = 0. Note that 
the period of cyclotron motion T and the cyclotron mass 
m c , given by Eq.Q, also have a singularity for the same 
reason here. 

Next we consider a|| y and a z-dependent modulation 
potential. The z-dependence of the modulation affects 
mainly the role of a real-space trajectory in the vicinity 
of the turning points related to kp^ x and — kp }X . Ac¬ 
cording to Eq. ([5]) the z-coordinate depends on k x and an 
electron sees different electric fields near these two points. 
As a consequence, the electric field averaged over the cy¬ 
clotron orbit is modified. The amplitude of oscillations 
does not drop to zero, as can be seen in the Fig{3] d), 
and the weights of frequencies /1 and / 2 in the frequency 
spectrum change. 

The main results of our model calculations are shown 
in the Figj4} The evolution of the Fermi contour calipers 
in the in-plane magnetic field B y is presented together 
with the values derived from the FFT of commensura¬ 
bility oscillations. The oscillations were calculated as¬ 
suming a z-dependent modulation potential that yielded 
the reduced electric field in the second layer. Here 
£{—d/ 2) = 0.4 x £{d/ 2) was chosen. 

Our calculation does not take into account the line 
broadening by the inaccuracy of instruments, by the fi¬ 
nite temperature and/or by electron scattering. On the 
other hand, the line broadening strongly influences inter¬ 
pretation of experimental data. Even in the most recent 
measurement^!] only one frequency was resolved in oscil¬ 
lation spectra measured both for a || x and a || y. For the 
case of a || x the good agreement with the frequency cor¬ 
responding to kp y max was found, for a || y the measured 
frequency is between /1 and / 2 mentioned above. 
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FIG. 4. The in-plane field dependence of the Fermi contour 
caliper dimensions. For a || x the kF, y max splits into kF, y min 
and ( kF, y ma x+kF, y min) /2 above the critical field B y for which 
the Fermi contour becomes concave. For a || y, the value 
kF,x max/2 appears abruptly at this critical field. An effective 
kF,x eff is derived from the weighted average of frequencies fi 
and fi corresponding to those vectors. 

In our calculation, which does not take into account 
the line broadening, the agreement between the Fermi 
contour calipers and the values derived from the FFT 
is very good for a || x. Therefore only the curves cor¬ 
responding to contour calipers are plotted in the Fig|4] 
as functions of B y . The kp, y max splits into k,F, y min 
and (kp, y max + ^F,y min)/2 above the critical field B y for 
which the Fermi contour becomes concave. The weight 
of the frequency related to ftp, y max dominates. 

The case of a || y is different. For B y below the critical 
value only the peak related to fi is observed in the FFT 
spectrum. The second peak corresponding to fi = /i/2 
appears with a very small weight at the critical field and 
dominates at high fields. To mimic the experimental re¬ 
sults, which resolve only one frequency, we plot an ef¬ 


fective kF,x eff, derived from the weighted average of fre¬ 
quencies fi and fi, except of the values based on the 
frequencies fi and f 2 themselves. 


V. SUMMARY AND CONCLUSIONS 

We have studied the commensurability oscillations in 
a quasi-two-dimensional electron gas modulated by a 
unidirectional periodic potential and subject to tilted 
magnetic fields with a strong in-plane component. A 
quasi-classical approach to the theory, that relates the 
magneto-resistance oscillations to the guiding-center drift 
was employed. A systematic discrepancy between the 
observed and the calculated elongation of the Fermi con¬ 
tour for the in-plane magnetic field parallel to the lattice 
vector, found in^, is explained by the in-plane magnetic- 
field distortion of the cyclotron orbits . 

Our arguments are as follows: First, the velocity of an 
electron moving along an orbit is no longer constant and 
depends on its position on the Fermi contour. Second, 
the cyclotron orbit is tilted, the electron moves from one 
side of a well to the other, and sees the different strengths 
of the one-side-modulation potential on the different well 
sides. 

The numerical simulation is based on the simple tight- 
binding model of a double well, the spin of electrons is 
not taken into account. 
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